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ABSTRACT 

We present our discovery of a new hybrid pulsating subdwarf B star, J08069+1527. The effec- 
tive temperature and surface gravity of 28,500±400K and 5.37±0.04dex, respectively, place 
this object inside the instability strip and also among other pulsating hot subdwarfs of a hybrid 
nature, right next to another fascinating star: Balloon 090 100001 . From this proximity, we an- 
ticipated this star could pulsate in both high and low frequency modes. Indeed, our analysis 
of photometric data confirmed our prediction. We detected two peaks in the high frequency 
region and two other peaks at low frequencies. In addition, the amplitude of the dominant 
mode is very high and comparable to the dominant peaks in other hybrid subdwarf B stars. 
Since this star is bright, we performed time-series low resolution spectroscopy. Despite a low 
signal-to-noise (S/N) ratio, we were able to detect the main peak from these data. All our 
results strongly indicate that J08069H-1527 is a high amplitude pulsating hot subdwarf B star 
of hybrid nature. By analogy to the other pulsating sdB star, we judge that the dominant mode 
we detected here has radial nature. Future stellar modeling should provide us with quite good 
constrains as p- and g-modes presented in this star are driven in different parts of its interior. 

Key words: stars: subdwarfs, asteroseismology. 



Hot subdwarf B (sdB) stars are horizontal branch stars with masses 
about 0.5 M0 and very thin, in mass, hydrogen envelopes. Their av- 
erage effective temperature and surface gravity are about 30000 K 
and log g ^ 5.5, respectively. Although the future evolution of the 
sdB stars to the white dwarf cooling track is generally accepted, 
their formation prior to the sdB stage is under debate resulting 
in s everal mechanisms that involve single-star or binary evolution 
fe.g. lD'CruzetaI.ll9 96': Han et al. 2002, 200?). 

The detection of pulsations in sdB stars opened a way to study 
their interiors and may help to understand their evolution prior to 
th e horizontal branch. F irst, short period oscillations were found 
bv lKilkennv et alj Il997l) in EC 14026-2647 (now officially named 
V391Hya). According to theoretical models, these pulsations are 
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attri buted to pressure mod es and are driven in t he outer part of the 
star ( ICharpinetet"ai1ll997h . Several years later, iGreen et alj j2003h 
announced another kind of variability in sdB stars. In this case, 
howev er, the pulsations wer e identified with gravity modes. Fol- 
lowing lpontaine et alj j2003l) they are driven deeper in the sdB stars 
than the pressure ones. Short-period sdBV (sdB Variable) stars are 
located at higher, while long-period sdBV stars at lower effective 
temperatures. 

The pulsation modes of the two kinds of sdBVs probe dif- 
ferent regions, and stars displaying both types are particularly in- 
teresting since their interiors can be better constrained. Such h 
brid pulsations were first found in DW Lyn (Schuhetal. 200j 
and another three objects with clear pulsations at both short and 
long periods have so far been rec ognized; B alloon090100001 
tearan et al.ll2005l: lOreiro et al.ll2005h. V391 Peg jLutz et al.ll2009l) 



and RAT 0455+1305 tearan & Fox-Machadd l2010^. Another star, 
CS 1246, has similar logTeft and \ogg and one high amplitude 
short periodicity detected. Unfortunately, because of noisy data, 
the detection threshold in the low frequency region is too high. 
The best data are in the r' filter, and the amplitude spectrum cal- 
culated from these data contains some excess signal in the low fre- 
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Table 1. Observational Log. 
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Figure 1. Light curve in a BG40 filter obtained in December 2009 at KPNO. 
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Figure 2. All data obtained on J08069+1527. The first night was taken at 
KPNO, while the second at Baker and the last two at Mercator. Note that 
data from the the first two nights were taken in a BG40 filter while during 
the last two we used a B filter (Geneva system). Each set of points stands 
for one night. Two vertical lines indicate coverage of spectroscopic data. 



quency region, but more data are needed to confirm if CS 1246 is 
a hybrid pulsator. This should encourage others to obtain more ob- 
servations of CS 1246. Extremely high precision photometry ob- 
tained with the Kepler spacecraft has revealed several long pe- 
riod sdBVs (Reed e t al. 20 10; Baran et al. 201 1) and one short pe- 
riod sdBV (Kawal eretalj|2010 ) that appear to also show hybrid 
modes, at least intermittently, that are so weak that they can only 
be detected using satellite data. An unusual hybrid pulsator in an 
eclipsing binary, 2M1938-I-4603, was found to show an exception- 
ally large number of weak (<0.5 mma) pulsation modes spanning 
all frequencies from t he long period to the short period regime 
testensen eTal]|2010j) . 

Typical amplitudes of sdBV stars showing high frequency os- 
cillations are below 1%, which corresponds to 10 mma. Most hy- 
brid stars have unusually high amplitudes (few per cent), mak- 
ing them extremely desirable. Since sdBV stars tend to be faint 
(12 > my > 17 mag) high amplitude pulsators are easier to ob- 
serve, even with 1 m cl ass telescopes. 

As was shown by iBaran & Fox-Machadol feOlO "). for hybrid 
stars the dominant modes in both frequency domains appear at very 
similar frequencies. These are ~ 0.3 mHz for gravity modes and ~ 
2.8 mHz for pressure modes. As for gravity modes, it may not be 
so surprising as the peaks are detected in a very small range, but for 
the pressure modes it is more astonishing. This may indicate that 
these stars are alike in internal structure and/or evolutionary status 
on the horizontal branch. In this paper we present our discovery of 
a new hybrid sdBV star with a high amplitude short period mode. 
We obtained data on four nights, including the discovery one. 



2 J080656.7-I-1527 

The target J08069-I-1527, has UV magnitudes FUV= 13.68, 
NUV = 13.86 in the GALEX All-sky Imaging Survey jMartin et all 
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I2OO5L AIS) observation obtained on February 2, 2006. After com- 
paring the UV magnitudes with visual photometry in the NO- 
MAD survey (5=13.88, 1^=14.23, J? =14.83), we flagged the 
target as a likely sdB star, and it entered a list of bright spec- 
troscopy targets that served as a poor weather backup for a pho- 
tometric observation run in the NOT s earch for pulsating sdBs 
in January 2009 (Ostense n et al.ll2010bh . The extracted spectrum 
was fitted to a model atmosphere spectrum, using the same pro- 
cedure as for the other sdBs observed during the search program. 
The derived effective temperature of 28,877(±202) K and surface 
gravity of 5.34(±0.03) dex along with its encouraging brightness 
made it a priority target for the n ext observing season. Recently, 
IVermes. Kawka & NemethI j201Cl ) have published their first results 
from a survey of UV-excess stars in the GALEX archive, and 
J08069-F1527 is included in their sample. Their determination of 
Toff and log g is fully consistent with ours, within the errors. 



3 TIME SERIES PHOTOMETRY 
3.1 Discovery data 

We performed time-series photometry of J08069-I-1527 on 30 De- 
cember 2009 using the KPNO 2.1m telescope. We used an Apogee 
CCD with a BG40 filter (wide band covering UBV filters range) 
with 12 sec exposure times. Seeing conditions allowed us to set 
2x2 binning resulting in very short readout times of about 1 sec. 
To determine if a star is variable or not, we used near-real-time 
software (developed at Missouri State University) during our ob- 
servations to monitor the light curve as it was obtained. Right after 
we started to collect data, it was obvious that this star changes its 
brightness and the amplitude of this variability is relatively high. 
In fact, the amplitude was so big that even from our real-time soft- 
ware it was easy to estimate the period and amplitude of the vari- 
ability. We estimated less than 10 minutes for the period and about 
20mmag for the amplitude. Based on its physical parameters, we 
assumed this variability to be caused by pulsation and realized that 
this is another (out of about seven known so far) high amplitude 
pulsating sdB star and decided to carry on observations on that 
night as long as possible. In total we collected 4 hours of data. 



3.2 Supplemental data 

After we found J08069-I-1527 to be an sdB pulsator with a high 
amplitude mode, we tried to obtain more time-series photometry to 
extend time coverage to lower noise level in the amplitude spec- 
trum and increase frequency resolution. Unluckily, another oppor- 
tunity to observe this star appeared only after 10 days. This long 
gap makes it harder to combine all data together since a window 
function gets too complicated and so does pre-whitening. Follow- 
up data were collected at Baker Observatory and with the Mercator 



© 2010 RAS, MNRAS 000,[T]l7] 



J08069+1527: A new hybrid sdBV. 3 



20 
15 
10 
_ 5 

E ° 
E 2 



0) 

1 



E 

< o 



1 1 1 1 1 — 


1 ' 1 ' 1 1 1 ' 

!i i 
1 1. ■ 


■| ' 

i jj 1 

I llilv S«, , , i 


1 


-1 f 


II M ' 

I 1 
,11 1 

II 1 1 
II 1 

1 II 1 

iUi u 


IL-I : , -^j^ 1^ , u 

-4f 


"ill 1 





2 3 4 5 

Frequency [mHz] 



Figure 3. Amplitude spectrum of KPNO data. The top panel shows the 
original data while the subsequent panels show pre-whitening as indicated 
by dashed lines. The solid horizontal line at 0.937 mma represents S/N = 4. 
The significant peak left in the amplitude spectrum appears at 2 hrs (half of 
the run length) and requires more data to be confirmed. 



telescope on La Palma. At Baker Observatory we used a Photo- 
metrics RS-1340 CCD with a BG40 filter. Exposure times were set 
to 30 sec along with 1 sec readout time. At the Mercator telescope 
we used the Merope CCD camera, which was recently upgraded 
with a new E2V frame t ransfer CCD with 2048 x 3074 illuminated 
pixels testensenll2oTo3) . Here we used a B filter (Geneva system) 
setting the data cadence period to about 31 sec. We include all this 
information in Tab[T] 

We calibrated all star images for instrumental effects (bias, 
dark, flat field) and extracted brightnesses of target and field stars. 
In the next step we calculated differential photometry using a 
few comparison st ars. All these t asks were done by means of the 
Daophot package ( IStetson|[l987l) with a graphical user interface 
(FinRed, developed by ASB). Because the field of view only cov- 
ers a few-arcminutes on the sky, the first order extinction (differ- 
ential extinction) was negligible and long term variations (on the 
scale of run duration) were de-trended by subtracting a cubic spline 
curve. We assume here that any change in brightness, on time scales 
of about 6 hours or longer, is not intrinsic to the star and is likely 
caused by non-perfectly removed extinction, either first or second 
order or other atmospheric effects. 

In the next step, treating all nights separately, we transformed 
all brightnesses to fluxes and subtracted a mean value. Resulting 
light curves were then centered around zero. Next, we analyzed 
these data by means of a Fourier analysis. We present our discovery 
light curve of J08069-I-1527 obtained at KPNO on 30 Dec 2009 in 
FigO One can note that besides short period variability some long 
term variability is also present in the light curve. We present all data 
in Fig|2] 
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Figure 4. Amplitude spectrum and the residual one calculated from data 
taken at Baker. The sohd horizontal hne at 2.84 mma represents S/N = 4. 



3.3 Amplitude spectrum 

Although we have data taken with three different telescopes, we 
decided not to combine them all together. It is because data from 
KPNO and Baker, taken in the same filter, are 10 days apart. When 
combining such data, we arrived with a higher noise level in the 
amplitude spectrum and a very complex window function. In the 
case of data from Baker and Mercator, although collected on con- 
secutive nights, they are not in the same filter. As a result, during the 
pre-whitening process, we can obtain residual (fake) peaks caused 
by different amplitudes and perhaps phases of modes observed in 
different filters. 

Since the telescope at KPNO is the biggest among those we 
used, data collected on that site are characterized by the smallest 
noise level in their amplitude spectrum. These data allow us to 
detect the highest number of peaks at S/N>4 used as our thresh- 
old. We detected two peaks in the high frequency region and two 
in the low frequency region. The dominant peak has a high am- 
plitude and a frequency similar to some o ther sdBV stars discov- 
ered so far. iBaran & Fox-Machadol ( |2010|) showed that all hybrid 
sdBV stars detected from ground have a main dominant high fre- 
quency peak at similar frequencies around 2.8 mHz. The main peak 
in J08069-F1527 is also located in the same region. What is more, 
its effective temperature and surface gravity is almost the same 
as other hybrids, particularly Balloon 090 100001. In addition, we 
found two peaks in the low frequency region. Precise estimation 
of their significance is rather hard since not all the signal can be 
pre-whitened from that frequency region. This increases the real 
noise level. A detection threshold used in our analysis is an aver- 
age through all frequencies in the amplitude spectrum, which might 
be a source of underestimation in the low frequency region. To 
estimate the noise level in that region more reliably, we used the 
adjoining 0.35-1 mHz region and assumed it to be valid for those 
two peaks. Although the detection threshold increased, they still re- 
mained significant. Amplitude spectra with all pre-whitening steps 
calculated from KPNO data are shown in Fig[3l 

In pre-whitening data taken at Baker, we could safely secure 
only one peak at the same frequency (within errors) as the domi- 
nant one from analysis of KPNO data. Others, if present are below 
the detection threshold. Only one peak detected here should not be 
surprising as the average noise level reaches 0.7 mma so only peaks 
higher than 2.8 mma could be detected with sufficient confidence. 
Looking at the result from KPNO, none but the dominant peak sat- 
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Figure 5. Amplitude spectram with pre-whitening steps calculated from 
data taken at Mercator. The solid horizontal at 1.65 mma line represents 
S/N = 4. 



isfies this condition. An amplitude spectrum as well as the residual 
one is shown in Fig|4] 

At Mercator we obtained data on two consecutive nights. Al- 
though here we have better resolution and more points collected, we 
could not detect more peaks compared to the KPNO data. There is 
some signal excess in the low frequency domain. Unfortunately, the 
noise level is too high to confirm any peaks in this region. We can 
only confirm two peaks in the high frequency region. The dominant 
peak we detected here has the same frequency as the one detected 
in KPNO and Baker data while the other has a similar frequency 
but barely exceeds the error limits. However, apart from amplitude 
errors, frequency errors are sensitive to the length of observations. 
That is why we may consider the frequencies derived from Merca- 
tor data as the more precise ones. On the other hand, after removing 
the dominant peak, the residual amplitude spectrum still contains 
a small peak at the frequency of the removed one (Fig|5). It can 
mean that this peak is physically unstable, or there may be a mode 
beating, or some instrumental effects that cause either frequency, 
or amplitude to be different on each night. From Tab|2]we can see 
that frequencies are not changing within the given eiTors. Compar- 
ing amplitudes is harder since at Mercator we used a different filter 
By comparing KPNO and Baker data, which were obtained in the 
same filter we can conclude that the amplitude has not changed over 
10 days (or is changing with periodicity of about 10 days). Trying to 
solve this issue we chopped Mercator data into nightly chunks and 
Fourier analyzed them separately. From the Fourier solution, we 
found that timing issues between nights is not bigger than 0.5 sec. 
However, the amplitudes from these two nights are substantially 
different. We derived 25.72(42) mma and 29.12(42) mma on the 
first and second nights, respectively. Since the residual of the pre- 
whitened data roughly equals the difference of the amplitudes from 
the two nights, we can conclude that the cause is the change in 
amplitude. 



4 TIME-SERIES SPECTROSCOPY 

We performed time-series spectroscopy on the night starting on 
7 Jan 2010. We used ALFOSC at the Nordic Optical Telescope 
(NOT) to obtain 472 low resolution spectra. We used grism#16 and, 
for all but the first few exploratory spectra, we used a 1.3 - arcsec 
wide vertical slit. Our dataset spans a total of 6 hours and 40 min- 
utes. 

We used a small window and 2x2 binning in the spatial di- 
rection to minimize readout-overheads; the window allowed for 
28 arcseconds of sky on either side of the stellar spectrum. To 
avoid under-sampling and to not smear the dominant variability too 
much, we set the exposure time to 30 sec, achieving a 42 - sec cycle 
time. In the course of our run we obtained spectra of a helium arc 
lamp to calibrate the wavelength. 

Unfortunately, since the brightness of this star is about 14 mag 
in the B filter only, we could not achieve very high S/N ratio. 
For single spectra we had typical S/N~15 in the 440-480 nm re- 
gion, while we reached S/N~150 in the average spectrum shown 
in Fig|6] 

We used standard IRAF procedures to reduce and calibrate the 
spectra, leading to 472 extracted spectra. We normalized the aver- 
age spectrum to unity, using a low-order polynomial fit. The result- 
ing spectrum was then fitted with a model spectrum, and the re- 
sulting model spectrum was used for guidance in the final normal- 
ization of the average spectrum. The stellar parameters we derived 
from the average spectrum (Teft, log g and log He in FiglTj were 
consistent for different normalization/rectification methods that ei- 
ther employ straight-line fits around each Balmer line or a multi- 
order polynomial fit that rectifies the complete spectrum at once. 
Then all individually extracted spectra were rectified using the fit 
that was needed to normalize the average spectrum. Finally the in- 
dividual spectra were scaled to get their continuum at unity, using 
a simple parabolic fit in the regions with no absorption lines. 

We derived Tcff =28,063±163 K, log 3 = 5.39±0.02 and 
log //e/// = -2.971 ±0.064 dex. These values are different than 
those derived from spectrum taken in January 2009 so we adopt 
an average value for all three parameters: Tcff = 28,500±400 K, 
logp = 5.37±0.05dex and logi^e/i? = -2.96±0.10. The reader 
should be aware, though, that systematic effects from the model 
grid can give shifts on the order of 2000 K and 0.2 dex in Tcff and 
log g, respectively. 

To look for radial velocity (RV) displacements we used the 
IRAF cross-correlation routine FXCOR, finding the radial-velocity 
shift of each individual spectrum with respect to that of the aver- 
age spectrum. We included only the wavelength regions around 7 
Balmer lines (H/3-HiO) in the cross-correlation. We fitted a Gaus- 
sian function to find the center of the cross-correlation function, 
which we adopted as a measure of the radial-velocity shift. As a 
result we obtained Heliocentric Julian Date along with RV shifts 
and their errors. We present the radial velocity curve in Fig [8] This 
result was subject to Fourier analysis in order to extract any period- 
icity from the data. The amplitude spectrum is shown in Fig|9] By 
means of non linear least square fitting, we detected only one peak 
which satisfies S/N>4. Its frequency, amplitude and phase with er- 
rors are: 2.8202(42) mHz, 6.1(1.0)km/s and 5.69(17) rad, respec- 
tively. The phase is given at HJD = 2455204.577483. We present 
the radial velocity curve folded with the period obtained in our anal- 
ysis in Fig[TO] The periodicity detected in a Fourier space barely 
meet our significance level as S/N = 4.55. 
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Figure 6. Average spectram of J08069+1527 obtained at NOT and calculated from 472 single low resolution spectra. Magenta line represents a model spectrum 
used for the final nomialization. 
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Table 2. Results of the pre-whitening process of all data treated separately by site. The phases are given at mean epoch (HJD) of each night: 2455195.960754, 
2455206.844458 and 2455209.065208 for KPNO, Baker and Mercator, respectively. The numbers in parentheses are the errors of the last digits. Gaps in the 
cells mean that no signal was detected either in Baker or Mercator data at frequencies detected in KPNO data. 
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Figure 8. Radial velocity displacement in function of time deiived from 
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Figure 10. Radial velocity displacement folded with the peiiod derived 
from our analysis. The phase goes over two cycles to better see detected 
periodicity. 

ani = mode, has a photometric amplitude in a B filter of 53 mma 
and a radial velocity (RV) amplitude of 19.2 km/s. The amplitude 
ratio of RV/B is 0.362. The dominant periodicity in J08069+1527 
has an amplitude of 27 mma in a B(Geneva) filter and RV = 6. Ikm/s 
which results in a ratio of 0.226. The difference in ratio is sub- 
stantial, however, the amplitude spectrum of J08069-I-1527 might 
be unresolved and the amplitude of the dominant mode affected by 
poor resolution. More data of J08069-I-1527 taken over several days 
may help to better compare results on these two interesting stars. 



5 SUMMARY AND CONCLUSION 

In this paper we presented our discovery of a new hybrid pulsating 
subdwarf B star, J08069-I-1527. Effective temperature and surface 
gravity place this object in the instability strip among other hy- 
brid subdwarf s. From our analysis of photometric data we detected 
four peaks, including two in the low frequency domain, confirm- 
ing our prediction. As the star is relatively faint, spectroscopic data 
is not sufficient to detect more than just the dominant mode we 
found in the discovery data. However, it proved that time-series 
spectroscopy at 2 m class telescope can support our analysis with 
radial velocity shifts down to 14 mag hot subdwarf stars, particu- 
larly when data are taken over several nights. 

As it turns out, physical parameters are not the only 
similarities J08069-I-1527 has to other hybrid sdBV stars. As 
shown in iBaran & Fox-Machadd ilOld in Fig. 9), some hybrids 
share a common frequency for their dominant variation and 
J08069-I-1527 also has this characteristic. It might be that these 
stars are in very similar evolutionary state and/or have very sim- 
ilar stellar structure. If so , we co uld use the result obtained by 
iBaran. Pigulski & O'Toolel ( |2008|) to constrain the degree of the 
dominant mode to 1 = 0. Balloon 090100001 has photometricall y 
constrained pulsation modes using multiplets jBaran et al.ll2009h . 
The dominant periodicity in Balloon 090 100001, associated with 
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